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Abstract 



We revisit the time evolution of a flat and non-flat direction system during in- 
flation. In order to take into account quantum noises in the analysis, we base on 
stochastic formalism and solve coupled Langevin equations numerically. We focus 
on a class of models in which tree- level Hubble- induced mass is not generated. Al- 
though the non-flat directions can block the growth of the flat direction's variance in 
principle, the blocking effects are suppressed by the effective masses of the non-flat 
directions. We find that the fate of the flat direction during inflation is determined 
by one-loop radiative corrections and non-renormalizable terms as usually consid- 
ered, if we remove the zero-point fluctuation from the noise terms. 



1 Introduction 



In the inflationary universe, quantum fluctuations of a scalar field are extended to super- 
horizon scale by the cosmic expansion. If the scalar field is massless, the variance of the 
superhorizon fluctuation grows as (H/2ir) 2 Ht, where H is the Hubble parameter during 
inflation and t is the cosmic time pQ. Inflaton, which is responsible for inflation, is one 
of the almost massless scalar field in the inflationary epoch and thus have quantum fluc- 
tuations. The slow-roll equation of motion for an inflaton is, however, usually described 
as a classical one. In order to analyze the effect of quantum fluctuations on the Inflaton 
dynamics, stochastic approach has been proposed [2j 13^1 . In the stochastic approach, one 
integrates all of the superhorizon modes, whose definition is somewhat artificial, to obtain 
so-called IR mode. The equation of motion for IR mode is described by coupled Langevin 
equations with noise terms. These noise terms represent the quantum "kicks" by horizon 
crossing modes and drive IR mode to evolve stochastically. 

In the context of supersymmetric standard models, there are many flat directions in 
scalar field space [5J. However, it is well-known that, in supergravity framework, an effec- 
tive mass of the order of Hubble scale is generically generated to a flat direction during 
inflation [6j [7]. Such a large effective mass can be avoided, if one assumes the D-term 
inflation scenario [6] or imposes a Heisenberg symmetry on Kahler potential [8]. Without 
tree-level Hubble-induced masses, one may have considered that one-loop radiative cor- 
rections and non-renormalizable terms determine the time evolution of a flat direction. 
Recently, Ref. [9] analyzed interacting systems which consist of flat and non-flat direc- 
tions without the tree-level Hubble-induced mass terms using the stochastic approach. 
In Ref. [9], it is insisted that the variance of the flat direction will saturate because of 
its effective mass in the potential arising from the couplings with the non-flat directions. 
However, in Ref. [9], the noise terms do not reflect the time evolutions of the effective 
masses. Thus, the conclusion in Ref. [9] can be essentially modified when we properly 
include the effective mass effects in the noise terms. In this study, we formulate the ef- 
fective mass effects in the noise terms and then analyze the time evolution of a flat and 
non-flat direction system. 

The construction of this paper is as following: in section [2j we review the formalism 
of stochastic approach to a single scalar field case. We also describe contributions of the 
zero-point fluctuation for the noise terms. In section [31 we formulate the noise terms 
including the effective mass of the scalar field. Then, in section HJ we consider a concrete 
system and show the results of the numerical calculation for coupled Langevin equations. 
Section |5j is devoted to conclusion. 

2 Stochastic approach 

In this section, we briefly review the stochastic approach to the time evolution of a scalar 
field [21 E]. In the following, we consider a real scalar field 4>(x) during inflation and take 

1 For applications to various inflation models, for example see Ref. [J| 
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Friedmann- Robertson- Walker metric: ds 2 = dt 2 — a 2 (t)dx 2 , where a(t) is the scale factor. 
For convenience, we also formulate contributions from the zero-point fluctuation for noise 
terms, which we use in the following analysis. 



2.1 Langevin equation 

Let us start with the equation of motion for a real scalar field <f)(x): 

where " ' " represents the time derivative d t and V(<f>) is a potential for <p(x). We define 
the mode expansion for <p(x) as 

/d 3 k 
(a k <p k (t) + a-k^(t)) e 4k ' x , (2) 

where the annihilation and creation operators a^, satisfy the following commutation 
relation: 

[a k ,4,] = ^(k-k'). (3) 

Here and hereafter, k is a comoving momentum. After linearizing the equation of motion 
Eq. (pQ), the mode function (pk(t) can be determined from the following equation: 

k 2 

(p k {t) + 3H<p k (t) + -jMt) + m 2 <p k (t) = 0, (4) 



a 2 



where the effective mass m is given by 



Here, (• • • ) is the spatial average inside the horizon during inflation. In the numerical 
calculation, we replace this spatial average by IR mode field value which is defined below. 
The solution for Eq. (j3J) is given by 



where Hi 1 \z) is a first kind Hankel function of order v and we have chosen the Bunch- 
Davis vacuum [10] . Here, v is given by v 2 = | — j^. We note that this mode function (Q 

satisfies (ph(t) — >■ -7=^e l ~^H in the short wave length limit -4? — > oo. 

Now, let us define IR mode for the scalar field <f>(x) and its conjugate tt(x) = 4>{x) as 

d 8 * 



/d k 
^ 3/2 6(ea(t)H - k) (a^ k (t) + a_ k y£(t)) e 

/d 3 A; 
j^-^6(ea(t)H - k) (a k cp k (t) + a_ k cp* k (t)) e 
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where 6{t) is the step function. Then, Eq. ([rj is reduced to the following coupled Langevin 
equations for IR modes [3J E]: 



$(x,t) = n(x,t) + s w (x,t), 

n(x,t) = -3ittl(x,t) - + s W( x ,t). 

Here the noise terms s^(x, i), s^(x, t) are given by 

d 3 k 



(8) 



s w (x,t) = ea(t)H 2 
s W (x,t) = ea(t)H 2 



(2tt) 3 / 2 

d 3 /c 
(2tt) 3 / 2 



5(ea(t)H - k) (a k ip k (t) + a_ k <^(i)) e 



8(ea(t)H - fc) (a k fe (t) + a_ k ^(t)) e 



ikx 



(9) 



:k-X 



Eq. (jHJ) is what we would like to derive in this subsection and is the basis of this study. 
Although we have considered a single real scalar field case, it is straightforward to derive 
Langevin equations for a multi- field system. In Section HI we analyze Langevin equations 
for a flat and non-flat direction system. 

In the numerical calculation, we need correlation functions between the noise terms, 
which are integrated over a short interval [t, t + dt] . The integrated correlation functions 
have the following forms: 



S w (r,t- dt) 



t+dt 



dtt 



t+dt 



dt 2 (vac| sW(x 1 ,ti)aW(x 2 ,t 2 ) |vac) 



H 



2tt 



£-) Hdtj (eaHr) ^e 3 |#«(e)| 2 



1 (v — real), 



(10) 



S w M;dt) 

/t+dt /-t+dt 
dti y dt 2 (vac| s (7r) (xi,ti)s w (x 2 ,t 2 ) |vac) 



rj2 \ 2 

—J Hdtj (eaHr) ^e 3 



3 

z/ 

2 



1 (z/ = real), 



(11) 



where r = |x x — x 2 | and jo(x) is a spherical Bessel function of order 0. Here, we have 
divided the expressions according to the cases where v is real and v — ifi (// > 0) is pure 
imaginary. 
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2.2 The zero-point fluctuation 

In this subsection, we formulate contributions from the zero-point fluctuation for the noise 
correlation functions. The mode function for the zero-point fluctuation is given by 

m™*°(t) = 1 _ p -ifdt'y/k*/a*+m* 

1 m 2 /H 2 1 ' ~ 2 — 1 (12) 

1+ ~- 



• zcro/jA 

Vk w = - 



2 (k/aH) 2 + m 2 /H 2 




Here, we have defined oo k {t) = ^Jk 2 + m 2 a 2 (t). Then, we obtain the following noise 
correlation functions: 



dN ^ 3 



S^ zero (r,N; dN) = —— Jo (eaHr) 



e 



S^ zero (r,N; dN) = Jo (eaHr) 



(2vr) 2JUV ' ^e 2 + m 2 /H 2 ' 
dN - 3 



1 Til 2 IH 2 \ 2 9 9,9 

l + » 2M ~ 2/m ) + e + m /H 
2e 2 + m 2 /H 2 J 



(2tt) 2JUV 7 y/e 2 + m?/H 2 

(13) 

Here, we have normalized the correlation functions by Hubble scale H as 5^(r, A/"; diV) = 
SW(r,t;dt)/iJ 2 , SW(r,N;dN) = S^(r, t; dt)/H 4 and N = Ht is the e-folding number 
(we need dimensionless quantities for the following numerical analysis). In particular, 
these expressions become simple for e 2 <C m 2 /H 2 : 



SW«*°( r ,N;dN) ~ -^ Jo (eaHr)e 3 ^ 
(2tc) 2 m 



Sl*y**°(r,N;dN) ~ (^o(eatfr)e 3 ^ (l + 
and for m 2 /# 2 = (i/ = 3/2): 



Am 2 /H 2 



(14) 



S(*y**°(r,N;dN) = £L Jo (eaHr)e 2 , 
(2n) 2 

S^°(r, N; dN) = -^L Jo (eaHr)e 2 (l + e 2 ) . 



(15) 



Now we mention the treatment of the contributions from the zero-point fluctuation. 
Although, as far as we know, we have no guiding principle to remove or to keep the 
zero-point fluctuation in the noise terms, we remove the zero-point fluctuation basically. 
The case where the zero-point fluctuation is left in the noise terms is discussed separately 
below. Fortunately, however, as we see in the next section, the contributions from the 
zero-point fluctuation are sufficiently suppressed by powers of e. Thus, as long as we take 
sufficiently small e, the effect of the zero-point fluctuation in noise terms are negligible. 
An exception is the case for large effective mass fh ^> H, where the zero-point fluctua- 
tion dominates the noise terms. In the next section, we discuss the effect of zero-point 
fluctuation on noise terms concretely. 
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3 Correlation functions between noise terms 



In this section, we formulate the noise correlation functions which we use for the numerical 
analysis in the next section. In particular, we take care of the effective mass dependences 
of the noise correlation functions. The noise correlation functions have quite different 

form between v = (-00 < §^ < |) and v = i/i = i\J%~ \ (§5 > f )• Thus, 

in this section, we formulate the noise correlation functions for each case. 



3.1 case : v is real 

Here, we consider the case where v - 
approximation formulae Eq. ( 13TT) . ( 132 
Hankel function for e<l: 



I - % is real (-00 < % < f). From the 
in Appendix, we obtain the following form of 



-oo< 



i 

~7I 



(lnf) : 



m 

w 



(16) 



Then, the dimensionless noise correlation function is given by 

S^(r,N;dN) - ^ - " 



dN 

(27T) 2 



j {eaHr) x 



2 2 T(v) 2 z-2v 
2vr t 



-OO < 



H 2 



< 2 



2eV+^(lnf) i (l-2 I 



o fh 2 ^ 9 
— H 2 — 4 



(17) 



To make our analysis easy, we naively extrapolate = 2 to = § in the second line 
in Eq. ( IT71) . We use this extrapolated formula for the numerical calculation in Section. HJ 

~ 2 

There is no special reason for choosing ^ = 2 as the boundary, however, we have checked 
that the numerical results in Section. H]is insensitive to ^ > 2 for e = 10~ 2 . Note that 
the zero-point contribution Eq. f|T3|) is not a leading order in e compared with Eq. f JTTj) . 
Thus, the treatment of the zero-point contribution does not matter in this case. 

In the same way, we approximate the noise correlation functions for the conjugate 
tt(x) (e <C 1) as 

am ^ /q \ 1 

S {lT \r,N;dN) 



dN n 3 

3o{eaHr) x - e 6 



(2vr) 2 
dN 



j (caHr) x 



3 

2 ~ 

2 2v r(u) 2 

1 — - — ( 

2tt 



v 1 F«(e 
(I 



7/ 



-00 < 



< 2 



2^ + S(lnf) 2 (l"2 



1/ 



9 m 2 9 
Z - tf 2 - 4 



(18) 



In the second line in Eq. (|T8|) for — 00 < ^ < 2, the expression becomes just for ^ = 
{y = |). With more accurate approximation, it should be j^^joieaHr) x • 



7T6" 2 

2 fc 
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-^^jo^eaHr) x e 4 . Since this is sufficiently small for e <C 1, we omit this e 4 contribution 
and use Eq. (ITS]) in the following numerical calculation. We note that the numerical 
results in Section. @]is insensitive to the boundary = 2 for e = 1CT 2 . Also, note that 
the zero-point contribution is not a leading order in powers of e compared with Eq. ffT8]) 
and thus the treatment of the zero-point contribution does not matter here. 



3.2 case : v = ifi, (fi > 0) 



Next, we consider the case where the effective mass is very large: v — i/j, — i<J — | 

(§2 > f )• In this case, from Eq. ( 135]) in Appendix, Hankel function for e 1 is approxi- 
mated as 

{ 2e 7r ^ ( coth tt[i _ Aire- 2w ^ p g^^j \ (fa 2 ^ 9\ 
' \ - »-'-":' ,Btr 'w) ^ a v ' (w) 

Then, the dimensionless noise correlation function is given by 

S^\r,N; diV) 

- — jo(eo^r) e"'" x ^ 3 \H^(e)\ 2 



(2tt) 2 juv y 2 

, , r f ^3 coth(7r/j)-l (fa 2 25 

diV . . rT , e 7^ I IT 2 " — T 



" X \ ^(coth( 27 r)-l) + ¥(lnf) 2 (l-f) (f<#<f). 

(20) 

Notice that the zero-point noise contribution given by Eq. (|14p is removed in the third line 
in Eq. (120]) . Even if we keep the zero-point contribution, the noise term S^(r,N;dN) 
(for e 1) is still small for large effective mass fh since S^(r, N; dN) is almost inversely 
proportional to fh. As in the case with real u, we extrapolate ^ = ^ to ^ = | in the 
second line in Eq. ( 120]) for the sake of easier analysis. We use this extrapolated formula 
in the numerical calculation in Section. HI We have checked that the numerical results in 
Section. @]is insensitive to this choice of the boundary = ^ as long as e <C 1. 

In the same way, we approximate the noise correlation functions for the conjugate 
tt(x) (e <C 1) as 



S w (r,iV; dN) 

= jo(eaHr)e- n ^ x - e 3 



(2vr 



2 



l-¥)H${e) + eH$_ l {e) 



dN f e 3 /i (l + fa) (coth(^) - 1) > f ) 

( 2 -) 2 X I ^S(ooth(27r) - 1) + g (In f ) 2 (1 - f) 



9 <r ™- 2 <r 25 

4 — H 2 — 4 



(21) 
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Here, the zero-point noise contribution given by Eq. (|14p is removed in the third line. 
In contrast to S ^ (r, N; dN) , the behavior of S^(r, N;dN) is quite different for large 
effective mass m if we keep the zero-point contribution. This is because S^\r, N; dN) oc 
e 3 m e~ 2m / H when we remove the zero-point contribution, while (r, N; dN) oc e 3 m 
when the zero-point contribution is included for m H. Notice that S^'(r, N; dN) and 
( r? ]V; dN) follow the ordinary uncertainty relation when the zero-point fluctuation 
contributions are included. In the numerical analysis below, we take e = 10~ 2 and thus 
the zero-point fluctuation contribution to the noise term never become important at least 
for iV < 10 3 . Now, we note that the numerical results in Section. H]is insensitive to the 
boundary g = f . The reason for the insensitivity to the above artificial boundaries 
= 2, 25 is that the noise terms with 2 < ^ < oo are so small (for e < 1) that these 
terms cannot dominate the behavior of IR mode. In the numerical calculation, we have 

- 2 

checked that the behavior of IR mode is determined by the noise terms only with <2 
fore = 10~ 2 . 



4 Numerical analysis for flat and non-flat direction 

systems 

In this section, we prepare for the numerical calculation of the Langevin equations. The 
model we consider here is taken from Ref. [9] as a concrete example. Then, we show the 
numerical results and discuss about the feature of the time evolutions of IR modes. 



4.1 Setup 

Here, we rescale Eq. ([8]) for the purpose of the numerical calculation. We rescale fields by 
Hubble scale H during inflation as 

$^$ = $/#, S W = s w /H 2 , 

n ->■ n = 11/ H 2 , s w s« = s« /H 3 , (22) 
V($) V($) = V(<S>)/H\ 

where bar shows the rescaled quantity. Using these dimensionless quantities, Eq. OH]) is 
rewritten as 

A$ (iV ) = n(iv) + #)(iv), 

= -saw + *> W , 

where we have omitted the x-dependence of the noise terms, since IR mode is almost 
spatially homogeneous inside the horizon. We note that the time argument is now con- 
verted to the e-folding number iV = Ht. Discretizing Eq. ( )23|) . we obtain the following 
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dimensionless Langevin equations: 



$(iV + diV) = $(iV) + n(iV)diV + f 



(<p) 



dv(HN))^ T , m (24) 



Tl(N + dN) = U(N) - 3Tl(N)dN d<S + ^ 

Here, the noise terms f^, £ w are the Gaussian random variables satisfying the variance 
given by Eqs. (TTTI) . ( IT8|) . (!20|) and (12T1) depending on the effective mass m. Although 
Eq. 021]) is derived for a single real scalar field case, generalization to the multi-field case 
is straightforward. In the next subsection, we show the numerical results for the coupled 
Langevin equations with three real scalar fields, where all noise correlation functions 
between different fields vanish. 



4.2 Numerical analysis 

Following Ref. [9], we consider the following scalar potential as a generic model for a flat 
direction and non-flat directions 

V = \\l (0 2 + h 2 ) e 2 + \g*K* + \gl (h* + 4e* - 4h 2 e 2 ) + ^, (25) 

where X e ,gi,g 2 are the coupling constants. Eq. (F25]) is motivated by the LH U direction <p 
in minimal supersymmetric standard model (MSSM). Here, the Hubble- induced effective 
masses [SI [7] are assumed to be absent. This happens within D-term inflation scenarios [6] 
or when imposing a Heisenberg symmetry on Kahler potential [8], for example. In Eq. ( T25|) . 
the last term is the non-renormalizable term and Mp ~ 2.4 x 10 18 GeV is the reduced 
Planck scale. For the potential Eq. (I25p . the Langevin equations are written as 

$(iV + diV) = <J(iV) + n(iV)diV + 

(tt2 _ \ 
\ 2 e E{N) 2 §{N) + 6— 2$(A^) 5 diV + ^\ 

E(N + dN) = E(N) + U e (N)dN + £ (e) , 
Tl e (N + dN) = f[ e (N) - 3U e (N)dN 

- ((\ 2 MN) 2 + (A 2 - g{) h{Nf) E(N) + 2g 2 1 E(N) 3 ) dN + £H 
h(N + dN) = h(N) + U h (N)dN + £ {h \ 

f[ h (N + dN) = S & (JV) - 3tl h (N)dN 

- ((A 2 - g 2 ,) E(N) 2 h(N) + i (gf + g 2 2 ) h(N)^ dN + 

(26) 

2 Do not confuse the "bar" symbol in e with in quantities $ etc. . e is the right-handed slepton and 
$ etc. are the dimensionless quantities. 
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Here, $ (E, h) and II (II e , 11^) are IR modes for the scalar field (e, h) and its conjugate. 
In the numerical analysis, we set all the fields at the origin initially: 



$(0) = E(0) = h(0) = 0, 

n(o) = n e (o) = i\(o) = o. 



(27) 



In the following, we omit the non-renormalizable term in Eq. f l2T)j) for a moment, which 
never becomes important within the scope of this study. 

The noise terms depend on the effective masses of the scalar fields. For the sake of 
convenience, we write down the effective masses: 



where rh^, rh e and fhh are the effective masses of the scalar fields 0, e and h, respectively. 
As we have mentioned below Eq. (jSJ), we have replaced the spatial averages by IR mode 
field values in Eq. ff28l) . In the numerical calculation, we read out the effective masses 
rh^,rh e and rhh, and then determine the noise terms step by step. 

In the numerical calculation, we take the time-step as dN = 10~ 2 . For generating 
random numbers, we utilize the Mersenne Twister method [11] . Then, we obtain Gaussian 
noises by Box-Muller method with the random numbers. The number of trials in our 
calculation is 5 x 10 4 . 

4.3 The results 

Here, we show the numerical results for Eq. ff26l) . In the numerical analysis, we take e = 
10~ 2 to guarantee the validity of approximation for Hanckel functions. We have checked 
that the result is almost independent of e < 1.0, if we remove the zero-point fluctuation 
from the noise terms. We also take the coupling constants as A e = g% — g% — 1.0 El- We 
note that, for these couplings, Eq. is reduced to 



In Fig. [TJ we show the time evolutions of the variances for IR modes. The red solid 
line, the green dashed line, and the blue dotted line show the IR modes for the direction 

3 Although these coupling constants A e = g± = gi = 1.0 are not realistic ones in MSSM, the generic 
feature of the flat direction IR mode can be studied. It is because, as we see in Fig. |3l the flat direction 
reaches an exactly flat direction at late time and coupling constants affect only to the relaxation time 
(and saturation value for h in Fig.[TJ. 



m\lE- = XlE(Nf, 

rhl/H 2 = A 2 $(iV) 2 + (A 2 - gl) h{Nf + 6^ 2 E(iV) 2 , 
™l/H> = (A 2 - g\) E(N) 2 + \ {gl + , 2 ) /.(iV) 2 , 



(28) 



ml/H* = XlE{N)\ 

m 2 e /H 2 = A 2 $(A0 2 + QglE(N) 2 , 



(29) 



m 2 h /H 2 = \(gl + gl) h(N) 2 . 
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Figure 1: The time evolutions of the variances for IR modes. The red solid line, the 
green dashed line, and the blue dotted line show the IR modes for the direction ("flat"), 
e ("nonflat 1"), and h ("nonflat 2"), respectively. Here, the coupling constants are A e = 
9i = 92 = 1.0 and e = 10~ 2 . 



("flat"), e ("nonflat 1"), and h ("nonflat 2"), respectively. Although all IR modes' 
variances are degenerate until the first one Hubble time (e-folds N — 1) pasts, they split 
completely after e-folds iV ~ 10 since each effective mass evolves differently. As we can 
see in Eq. fT26|) . non-flat direction h is decoupled from others for A e = g\. Thus, the 
time evolution of IR mode for h is determined by the self-coupling term h 3 and the noise 
terms. From this self-coupling, the effective mass (fnfy/H 2 = 3(h 2 ) <C 1 is generated. 
The saturated value for h (~ 0.11) in Fig. [JJ can be understood by the formula for the 
variance with rh <C H [1J: (h 2 ) = %^r~m ■ Using this formula, we obtain the averaged 

value (h 2 ) = y/TJ (8tt 2 ) ~ 0.11 which is consistent with the saturated value in Fig. [TJ On 
the other hand, the flat direction <p and the non-flat direction e couple with each other. 
Since e has a self-coupling term and </> does not, the effective mass for e become lager 
than the one for (p. Thus, although the variances of their IR modes are degenerate at 
first, they split after the self-coupling of e becomes non-negligible in the effective mass. 
At iV ~ 10 3 , e becomes too massive and cannot have a large fluctuation. In contrast, 
eventually becomes as an exactly flat direction. See also Fig. |3] for the comparison with 
the case when is exactly free. 

In Fig. El we show the averaged variances for the integrated noise correlation functions 
(SW) (left panel) and (right panel). The noises are integrated for intervals dN = 

10~ 2 . Although the noises (S^'} for and e are degenerate at first, (S^) for eventually 
approaches the exactly flat case and (S^) for e decreases rapidly. This is consistent 
with Fig. [TJ We note that (S^) has little effect on the numerical calculation, since 
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Figure 2: The averaged variances for the integrated noise correlation functions (S^) 
(left panel) and (S^) (right panel). The noises are integrated for intervals dN = 1CT 2 . 
The lines, coupling constants and e are the same as in Fig. CD 




e-folds N 

Figure 3: Comparison with flat direction IR modes for various cases. The red solid line 
is the same as in Fig. [T] and the green dashed line is for A e = 0, g\ = #2 — arbitrary. The 
blue dotted line is the same as Ref. but with e = 10~ 2 for X e = g 1 = g 2 = 1.0 where 
the massless noise terms are used (1/ — 3/2 for <p,e, h). 
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In Fig. El we show the time evolution of flat direction IR modes for various cases. The 
red solid line is the same as in Fig. [IJand the green dashed line is the exactly flat direction 
case (A e = 0, g\ = gi — arbitrary). The blue dotted line is the same as Ref. [9] but with 
e = 10 -2 for A e = gi — g% = 1.0 where the massless noise terms are used iy = 3/2 for 
0, e, h). From Fig. [3j we confirm that our result is much different from the one in Ref. [9]. 
This is because Ref. jH] does not include the effective mass effects in the noise terms. When 
the e-fold number is small, our result and Ref. [§] are consistent with each other since all 
of the effective masses are sufficiently small. However, after the effective mass for non-flat 
direction e becomes large N ~ 10, the flat direction becomes more flat and e becomes 
more non-flat in our formulation. Then, the flat direction eventually approaches to an 
exactly flat direction in our case. Thus, our conclusion is different from Ref. [2] : even if 
some non-flat directions prevent a flat direction to go away from the origin, the flatness 
for the flat direction eventually recovers and the variance of flat direction increases as 
large as the exactly flat case. 

Now, we comment on the relevance of one- loop radiative corrections [8j [12] . Since the 
tree-level effective mass of the flat direction is highly suppressed at last, one-loop radiative 
correction \rh^\ ~ 0.01A 2 if 2 suggested in Refs. [HJ [T2] generally takes over the tree-level 
effective mass. Thus, it is reasonable to say that the time evolution of a flat direction 
is determined by one-loop radiative corrections and non-renormalizable terms. This is 
the same situation as one may usually have considered. We note that the zero-point 
fluctuation in the noise terms has been removed in the above argument. The case where 
we include the zero-point fluctuation is discussed in the next subsection. 

Finally, in Fig. HI we show the time evolutions of the probability density functions 
for the flat direction and the non-flat directions e, h. We construct the probability 
density functions by dividing the realizations into bins. We confirm that the flat direction 
diffuses as time goes. As the effective mass for approaches to for iV > 10, 
approaches to the exactly flat direction and will continue to diffuse for iV > 10 3 . On the 
other hand, although the non-flat direction e diffuses at first, e concentrates on the origin 
after iV = 10 (the values in Fig. H]for e at N = 10 2 , 10 3 are about 20,40, respectively). 
This is because e gets a large effective mass, namely, the noise term for e is small and the 
potential enforces e to go back to the origin for N > 10. This is consistent with Figs. [1] 
and |5J The probability density function for h does not change for N > 10. This is because 
the variance of h saturate at N ~ 10. Since the potential force is zero at the origin, the 
stable point is driven to non-zero h by the noise terms. This feature is also seen in the 
probability density function for e in Fig. H]at N = 10, 10 2 and 10 3 . 

4.4 Discussions 

Now, let us comment on the zero-point fluctuation contribution. If we include the zero- 
point contribution in the noise terms in Eq. f l2T|) . the noise term 5 , ^ 7r - ) (r, N; dN) oc e 3 m 
(for fh ^> H) eventually becomes large and induces a non-negligible effective mass to 
the flat direction through the coupling A 2 $ 2 £" 2 after many Hubble time past. Then, the 
variance of the flat direction will saturate at last. This means that the time evolution 
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Figure 4: The time evolutions of the probability density functions. The top-left, top- 
right, bottom- left and bottom- right panels are for flat direction 0, flat direction (log 
scale), non-flat directions e and h, respectively. The red solid line, the green dashed line, 
the blue dotted line , the purple thin dotted line and the emerald dash-dotted line are the 
probability density functions for N = 0.1,1, 10, 10 2 and 10 3 , respectively. The coupling 
constants and e are the same as in Fig. [TJ 
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of the flat direction will be quite different from the case where we remove the zero-point 
contribution. Fortunately, however, the existence of one-loop radiative corrections make 
the fate of flat direction during inflation not so different from the one without the zero- 
point contribution, as long as the renormalization scale is chosen appropriately. Anyway, 
we do not know which is the correct way to treat the zero-point fluctuation in the noise 
terms. 

5 Conclusion 

In this study, we have analyzed the time evolution of a flat and non-flat direction system 
governed by coupled Langevin equations during inflation. We have taken into account the 
effective mass effects on the noise terms. In the analysis, we have removed the zero-point 
fluctuation contributions from the noise terms. As the flat direction goes away from the 
origin, the effective masses of non-flat directions coupled to the flat direction become 
large. As a consequence, such a non-flat direction cannot have large fluctuations and 
cannot block the growth of the variance for the flat direction. Thus, the tree-level effective 
mass of the flat direction is eventually highly suppressed. The time evolution of the flat 
direction is, then, determined by one-loop radiative corrections and non-renormalizable 
terms as one have usually considered. We have also discussed the case where the zero- 
point contributions are included in the noise terms. In this case, the variance of the flat 
direction will saturate at last in the tree- level argument. However, little is known how 
to treat the zero-point fluctuation in the noise terms, while it is important for a massive 
field with the effective mass ^> H . 
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Appendix: Approximation forms for Hankel functions 

In this appendix, we describe the approximation formulae for Hankel functions H„ (e) (e <C 



1), which we use in this study. When v = J | — ^ is real (— oo < ^2 < §), we can use 
the following well-known approximation formula for Hankel function |13j : 





(30) 
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where e < 1 and T(u) is Gamma function. Its square is given by 



\H?\e)\^ 2 -^le-^. (31) 

We cannot use Eqs. fl3"U|) . ( IBTj) for z/ ~ and v = ifi. When v = 0, Hankel function is 
given by 

iT^(e) ~ 1 -i|ln| (e«l), 
|^ 1} (e)| 2 ^ — (in-) (e«l). 

On the other hand, when v = i/i (// > 0), we can apply the following approximation by 
using the expression for Bessel function with complex v [I] 

-1) 

V(k + 1YT 

fc=0 



" u £-< r(fc + i)ri> + * + 1) I2 



(33) 

Using this expression, we find the following approximation formula for Hankel function: 
HP(e) = J u (e) + iNM(e) 

= T ^Me)-e-J-Ae)) (34) 
2 ( 1 e--i^fir1 (e«l). 



l_ e 2mu \Y{l + v) \2J r(l-z/)V2 
Then, we find for v = %\i (fj, > 0): 

... A/1 p -2tt/i p 2iMlnf 



i-e- 2 ^| 2 V|r(i + i/x)| 2 |r(i-^)| 2 r(i + z/i) 2 

/cothvr/i 4vre~ 27rM 



Re- 



3 2i/x In ; 



(35) 



7T V A* (l-e~ 2 ^) 2 r(l + v/) 2 

References 

[1] A. D. Linde, Phys. Lett. B 116, 335 (1982): A. Vilenkin and L. H. Ford, Phys. Rev. 
D 26, 1231 (1982): A. A. Starobinsky, Phys. Lett. B 117, 175 (1982). 

[2] A. A. Starobinsky, in "De Vega, H.j. ( Ed.), Sanchez, N. ( Ed.): Field Theory, Quantum 
Gravity and Strings", 107-126. 

[3] M. Sasaki, Y. Nambu and K. -i. Nakao, Nucl. Phys. B 308, 868 (1988). 



15 



[4] K. -i. Nakao, Y. Nambu and M. Sasaki, Prog. Theor. Phys. 80, 1041 (1988); Y. Nambu 
and M. Sasaki, Phys. Lett. B 219, 240 (1989); S. Mollerach, S. Matarrese, A. Ortolan 
and F. Lucchin, Phys. Rev. D 44, 1670 (1991); J. Martin and V. Vennin, Phys. Rev. 
D 85, 043525 (2012) jarXiv:1110.2070l [astro-ph.CO]]. 

[5] T . Gherghetta, C. F. Kolda and S. P. Martin, Nucl. Phys. B 468, 37 (1996) 
[hep-ph/9510370| . 

[6] E. D. Stewart, Phys. Rev. D 51, 6847 (1995) |hep-ph/9405389| . 

[7] M. Dine, L. R andall and S. D. Thomas, P hys. Rev. Lett. 75, 398 (1995) 
|hep-ph/9503303l ; Nucl. Phys. B 458, 291 (1996) | hep-ph/9507453| . 

[8] M. K. Gaillard , H. Murayama and K. A. Olive, Phys. Lett. B 355, 71 (1995) 
[hep-ph/9504307| . 

[9] K. Enqvist, D. G. Figueroa and G. Rigopoulos, JCAP 1201, 053 (2012) 
jarXiv: 1 109. 3021 [astro-ph. GO]]. 

[10] T. S. Bunch and P. C. W. Davies, Proc. Roy. Soc. Lond. A 360, 117 (1978). 



[11] Mersenne Twister, http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/mt.html 



[12] B. Garbrecht, Nucl. Phys. B 784, 118 (2007) |hep-ph/0612011|. 

[13] N. N. Lebedev, R R. Silverman (translator), "Special Functions and Their Applica- 
tions", Dover Publications (1972). 

[14] K. Enqvist, K. W. Ng and K. A. Olive, Nucl. Phys. B 303, 713 (1988). 



16 



